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ABSTRACT 



Context. Open clusters older than ~ 4 Gyr are rare in the Galaxy. Affected by a series of mass-decreasing processes, 
the stellar content of most open clusters dissolves into the field in a time-scale shorter than ~ 1 Gyr. In this sense, 
improving the statistics of old objects may provide constraints for a better understanding of the dynamical dissolution 
of open clusters. 

Aims. Our main purpose is to investigate the nature of the Globular cluster candidate FSR 1716, located ai £ — 329.8° 
and b — —1.6°. We also derive parameters of the anti-centre open cluster Czernik23 (FSR 834). Both objects have been 
detected as stellar overdensities in the Froebrich, Scholz & Raftery star cluster candidate catalogue. 
Methods. The analyses are based on near-infrared colour-magnitude diagrams and stellar radial density profiles. The 
intrinsic colour-magnitude diagram morphology is enhanced by a field-star decontamination algorithm applied to the 
2MASS J, H, and Ka photometry. 

Results. Isochrone fits indicate that FSR 1716 is more probably an old (~ 7 Gyr) and absorbed (Av = 6.3 ± 0.2) 
open cluster, located ~ 0.6 kpc inside the Solar circle in a contaminated central field. However, we cannot rule out 
the possibility of a low-mass, loose globular cluster. Czernik 23 is shown to be an almost absorption-free open cluster, 
~ 5 Gyr old, located about 2.5 kpc towards the anti-centre. In both cases. Solar and sub-Solar {[Fe/H] ~ —0.5) 
metallicity isochrones represent equally well the stellar sequences. Both star clusters have a low mass content 200 Mq) 
presently stored in stars. Their relatively small core and cluster radii are comparable to those of other open clusters 
of similar age. These structural parameters are probably consequence of the several Gyrs of mass loss due to stellar 
evolution, tidal interactions with the disk (and bulge in the case of FSR 1716), and possibly giant molecular clouds. 
Conclusions. Czernik 23, and especially FSR 1716, are rare examples of extreme dynamical survivors. The identification 
of both as such represents an increase of ~ 10% to the known population of open clusters older than ~ 4 Gyr in the 
Galaxy. 
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1. Introduction 

The Galaxy is an aggressive environment to star clusters in 
general, the open clusters (OCs) in particular. These stellar 
systems are continually harassed by a series of dynamical 
processes such as mass loss associated to stellar evolution, 
mass segregation and evaporation, tidal interactions with 
the Galactic disk and bulge, and collisions with giant molec- 
ular clouds. Combined over time, such processes tend to 
accelerate the dynamical evolution, which produces signifi- 
cant changes in the cluster structure and eroded mass func- 
tions. Eventually, most OCs end up completely dissolved in 
the Galactic stellar fi eld or as poorly-populated remnants 

(|Pavani fc Bical[2007l and references therein) . 

Theoretical (e.g. ISpitzeil 119581: 
Lamers fc Giele sI 120061) . N-body (e.g . 



~ 1 Gyr. Around this region, t he disruption t ime-scale 
depends on mass as t^is ^ Af°-^^ (jLamers fc GieleSi2006i) , 
which for clusters with mass in the range 10^ — lO^M© 
corresponds to 75 ^ idis(Myr) ^ 300. In general, the 
effect of the relevant dynamical processes is stronger 
for the O Cs more centra l ly loca ted and the low-mass 
ones (see iBonatto fc Bical l2007al for a detailed discus- 
sion on disruption efi^ects and time-scales). Indeed, OCs 
older than ~ 1 Gyr are preferentially foun d near the 
So lar circle and in the outer Galaxy (e.g. iFriell 119951 : 
iBonatto et al.|[2006af ). where the frequency of potentially 
damaging dynamical interactio ns with giant molecular 
cloud s and the disk is lower (e.g. Salaris. Weiss fc Percivall 



Baumgardt fc Makinol |2003| iG oodwin fc BastianI 120061: 
Khalisi. Amaro-Seoan c & SpurzcirJ_ 20071). a nd observa- 



tiona l (e.g. van den Bergh il957. : lOort 119581 : Ivon Hoerneil 
1X9581 : iPiskunov et al.l 2007 ) evidence indicate that the 
disruption-time scale near the Solar circle is shorter than 
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120041: lUpgren. Mesrobian fc Kerridgd 119721 ). Disruption 
efficiency increases critically towards the Galactic cen- 
tre, to the point that the inner {Rgc ^ 150 pc) tidal 
fields can dissolve a massi ve star cluster in ~ 50Myr 
(jPortegies Zwart et al.ll2002[ ). 

The above aspects considered, the natural expectation 
is that only a small fraction of the OCs can reach old ages, 
and that the successful ones should be preferentially found 
at large Galactocentric distances. In fact, of the « 1000 OCs 
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with known age listed in the WEBDAQ database, 180 are 
older than 1 Gyr, and only 18 2 %) are ol der tha n 4 Gyr 
(see also lOrtolani et alj 120051 : IOrtolani~ Bica fc Barbuvl 
[20051) . Not surprisingly, most of the OCs older than 1 Gyr so 
far identified are located outside the Solar circle (see, e.g., 
the spatial distribution of OCs of different ages in Fig. 1 of 
iBonatto fc Bicall2007a|) . 

In this context, it is naturally expected that the discov- 
ery and derivation of astrophysical parameters of old OCs 
will better define the OC parameter space. Thus, a better 
understanding on the dynamical survival rate of star clus- 
ters in the tidal field of the Galaxy can be reached. Such 
parameters, in turn, can be used in studies of star forma- 
tion and evolution processes, dynamics of N-body systems, 
and the geometry of the Galaxy, among others. 

In the present paper we study two stellar overden- 
sities listed in the star cl u ster candidate catalogue of 
iFroebrich. Scholz fc Raftervl (|2007D . which turn out to be 
very old star clusters. They are FSR 1716 and FSR 834. 
The present work employs near-IR J, H , and Kg photom- 
etry obtained from the 2MAS^ Point Source Catalogue 
(PSC). The spatial and photometric uniformity of 2MASS, 
which allow extraction of wide surrounding fields that pro- 
vide high star-count statistics, are important to derive clus- 
ter parameters a nd probe t he nature of stellar overdensities 
(e.g. Bica, Bona tto fc Cam argo 2008). For this purpose we 
have developed quantitative tools to statistically disentan- 
gle cluster evolutionary sequences from field stars in colour- 
magnitude diagrams (CMDs). Basically, we apply (i) field- 
star decontamination to quantify the statistical significance 
of the CMD morphology, which is important to derive red- 
dening, age, and distance from the Sun, and (ii) colour- 
magnitude filters, which are essential for intrinsic stellar 
radial density profiles (RDPs), as well as luminosity and 
mass functions (MFs). In particular, field-star decontami- 
nation constrains more the age and dist ance, especially for 
low-latitude OCs (jBonatto et al.ll2006af) . 

This paper is organised as follows. Sect, [^contains ba- 
sic properties and reviews literature data (when available) 
on both star cluster candidates. In Sect. [3] we present the 
2MASS photometry and build the stellar surface-density 
distribution in the direction of both objects. In Sect. [3] 
we build CMDs, discuss the field-star decontamination al- 
gorithm, and provide tests to the old age of both clus- 
ters. In Sect. [5] we derive cluster fundamental parameters. 
Sect.[6]describes cluster structure by means of stellar RDPs. 
In Sect. [7] mass functions are built and cluster masses 
are estimated. In Sect. E] aspects related to the structure 
and dynamical state of the present clusters are considered. 
Concluding remarks are given in Sect. [HI 

2. FSR 1716 and FSR 834 as stellar overdensities 



The catalogue built by IFroebrich. Scholz fc Raftervl ()2007D 
includes 1021 star cluster candidates (hereafter FSR ob- 
jects) for Galactic latitudes |6| < 20° and all longitudes. 
The targets were selected by an automated algorithm that 
basically identifies small-scale regions that present stellar 
overdensities, applied to the 2MASS database. The over- 



^ obswww.univie.ac.at/webda - \MermiUiod & Paunzeri (20031) 



^ The Two 
data release 



Micron All Sky su 

Skrutski e et al.1 Il997f ). 



Survey, All Sky 
available at 



densities are classified according to a quality flag, '0' and 
'1' for the most probable star clusters, while '5' and '6' 
may be related to field fluctuations. Based on a combina- 
tion of the number of cluster stars (corrected to a common 
mag nitude limit), the core radius and t he central star den- 
sity, IFroebrich. Scholz fc Raftervl (|2007[ ) could discriminate 
known globular clusters (GCs) from OCs. With this crite- 
rion applied to the overdensity catalogue, they found 1012 
OC and 9 GC candidates. 

Several studies have explored the FSR catalogue 
with different approaches, with results that reflect 
the importance of such catalogue. The recently dis- 
cover ed GCs FSR 1735 (IFroebrich. Meusinger fc Scholz] 
I2007D and FSR 1767 (iBonatto et all 120071) . and the 
probable GCs FSR 584 (iBica et alJ 120071) and FSR 190 
dFroeb rich. Meusinger fc DavisI 12007 ). are clear examples 
of the fundamental role played by the FSR catalogue to 
improve the statistics of very old star clusters. Indeed, 
FSR 1735 and FSR 1767 are the most recent additions 
to the Galactic GC population, a number that presently 
amounts to ~ 160 members (e.g. Bonatto & Bica 2008al). 
On the other hand, the decontamination algorithm (Sect. 2]) 
applied to the 2MASS photom etry of some FSR GC c andi- 
dates, has shown that FSR 89 d^natto fc Bical l2007a|) and 
FSR 1603 (|Bica fc Bonattol I2008D are open clusters of age 
~ 1 Gyr, while FSR 17 55 appears to be an embedded OC 
(|Bica fc Bonattdl2008l ). 

During the course of a close investigation of the over- 
densities, we noted that the CMDs of FSR 1716 and 
FSR 834 show features typical of old stellar systems. 
Indeed, FSR 1716 was clas si fied a s a GC candidate by 
IFroebrich. Scholz fc RafteT^ (|2007l ). FSR 1716, at i = 
329.79° and b = —1.59° is projected not far from the 
Galactic centre, which implies the presence of significant 
contamination by bulge stars. FSR 834, on the other hand, 
is a disk object projected towards the anti-centre {£ = 
180.55° and b = -^0.82°). It has the poorly-stu died OC 
Czern ik 23 (hereafter Cz 23) as optical counterpart. ICzernild 
(|l966l ) measured a diameter of 5' and esti mated tha t Cz 23 
contains about 40 member stars, and lRuprech"3 (jl966| ) 
which classified it as III Ip. 

Figure [1] (left panel) shows a 4' x 4' 2MASS A% im- 
age of FSR 1716, where a significant concentration of stars 
is superimposed on a relatively crowded field, as expected 
from such a central direction. For Cz 23 we show in the 
right panel a 7' x 7' XDS^ R band image. In this case the 
cluster can be easily seen against a relatively clean field. 

Table [T] provides fundamental data on both objects. 
Coordinates from Frocbrich, Scholz fc Raftcry (2003) are 
given in cols. 2 to 5; their quality flag in is col. 6, and the 
core and tidal radii measured in the H band are given in 
cols. 7 and 8. The age, central reddening, distance from 
the Sun, Galactocentric distance, and the components pro- 
jected onto the Galactic plane derived in the present study 
(Sect. [5]) are given in Cols. 9 to 16. 



3. 2MASS photometry 

2MASS photometry in the J, H, and Kg bands was ex- 
tracted in circular flelds of extraction radius i?cxt centred 
on the coordinates of the objects (Table [T]) by means of 



http :/ / www. ipac. caltech.edu/ 2mass /releases /al lsky/\ 



^ Extracte d from the Canadian Ast ronomy Data Centre 
(CADC), at http://cadcwww.dao.nrc.ca/ 
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Fig. 1. Left panel: 4' x 4' 2MASS if^ image of FSR 1716. Image provided by the 2MASS Image Service. The small circle 
indicates the central coordinates (cols. 2 and 3 of Table[T]). Right panel: 7' x 7' XDSS R image of Cz 23. Figure orientation: 
North to the top and East to the left. 



Table 1. General data on the clusters 









FSR2007 














This paper 






Cluster 


a(2000) 


(5(2000) 


£ b 


Q 


Rc.u 


Rt,H 


Age 


Av 


dQ 


Rgc 




Ygc 


Zgc 




(hms) 


(°"') 


n n 




{') 


(') 


(Gyr) 


(mag) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(kpc) 


(1) 


(2) 


(3) 


(4) (5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


FSR1716t 


16:10:33 


-53:44:12 


329.79 -1.59 


2 


1.2 


5.9 


7.0 + 1.0 


6.3 + 0.2 0.8 + 0.1 


6.6 + 0.1 


-6.6 + 0.1 


-0.37 + 0.04 


-0.02 + 0.01 


FSR1716t 


16:10:33 


-53:44:12 


329.79 -1.59 


2 


1.2 


5.9 


12.0 + 2.0 


6.3 + 0.4 


2.3 + 0.3 


5.4 + 0.3 


-5.3 + 0.3 


-1.13 + 0.17 


-0.06 + 0.01 


Cz 23 


05:50:07 +28:53:28 


180.55 +0.82 


1 


0.8 


4.2 


5.0 + 1.0 


0.0 + 0.1 


2.5 + 0.1 


9.7 + 0.2 


-9.7 + 0.1 


-0.02 + 0.01 


+0.04 + 0.01 



Ta ble Notes. C oordinates (cols. 2 to 5 ), quality flag (col. 6), and core and tidal radii (cols. 7 and 8) measured in the H band are from 
IFroebrich. S cholz & Rafter^ I I2007D : Col. 10: reddening tow ards the cluster's central region (Sect. [Sjl. Col. 11: distance from the Sun. 
Col. 12: cluster Galactocentric distance for _Rq = 7.2 kpc lIBica et al.ll200^ ). Cols. 13-15: coordinate components projected onto the 
Galactic plane. Cz23 is the optical counterpart of FSR 834. Parameters of FSR 1716 are derived for the OC (f) or globular cluster (|) 
interpretation (Sect. [Sjl. 



VizieRS. Previous studies with OCs in different environ- 
ments (Sect.[T]) revealed that in the absence of a neighbour- 
ing populous cluster and important differential absorption, 
wide extraction areas provide suitable statistics for a consis- 
tent characterisation of the field stars in terms of colour and 
magnitude. Based on this premise we adopted an extrac- 
tion radius of i?cxt = 30', which is beyond the cluster radius 
(Sect, inland col. 9 of Table[3]) of the present objects. For de- 
contamination purposes, comparison fields were extracted 
within wide rings located beyond the cluster radii. As a pho- 
tometric quality constraint, the 2MASS extractions were 
restricted to stars (i) brighter than those of the 99.9% 
Point Source Catalogue completeness limit[3 in the clus- 
ter direction, and (ii) with errors in J, H, and Ks smaller 
than 0.3 mag. The 99.9% completeness limits refer to field 
stars, and depend on Galactic coordinates. In the present 
cases, the fraction of stars with J, H, and Kg uncertainties 
smaller than 0.06 mag is « 80%. A typical distribution of 

* fhttp: //vizier. u-strasbg.fr/viz-bin/VizieR ?-source=II/246 
^ \http://www.ipac.caltech. edu / 2mass /releases /allsky/ doc/ 



uncertainties as a function of magnitude, for objects pro- 
je cted towards the central parts of the Galaxy, can be found 
in lBonatto &: Bical (|2007b[ ). Reddening transformations use 
the relations Aj/Ay = 0.276, Ah I Ay = 0.176, AxJAy = 
0.118 , and Aj = 2.76 x E{J - H) (Dutra. Santi ago fc Bical 
|2002[ ). for a constant total-to-selective absorption ratio 
Ry = 3 .1. These ratios were derived from the extinction 
curve o f ICardelh. Clavton fc Mathii JTMl)- 

As a first step to understand the nature of both over- 
densities, we show in Fig. [2] the spatial distribution of the 
stellar surface-density as measured with 2MASS photom- 
etry. In both cases we compute the surface density (cr, in 
units of stars arcmin"^) in a rectangular mesh with cells 
of dimensions 4' x 4'. The mesh reaches total offsets of 
Aa| = |Ai5| ~ 20' with respect to the centre, in right as- 
cension and declination. Since the cluster radius of both ob- 
jects is <i 6' (Table [3]), most of the cluster structure should 
be contained in the central cell. Because of the important 
contamination by field stars, the surface-density built with 
the observed (raw) photometry of the centrally projected 
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FSR1 71 6 - Raw Cz23 (FSR0834) - Raw 




FSR1 71 6 - Clean Cz23 (FSR0834) - Clean 




Fig. 2. Stellar surface-density cr(stars arcmin ^) of 
FSR 1716 (left panels) and Cz23 (right). The curves were 
computed for a mesh size of 4' x 4', centred on the coordi- 
nates in Table [TJ The observed (raw) and field-star decon- 
taminated photometry are shown in the top and bottom 
panels, respectively. 



OC FSR 1716 is very irregular (top- left panel), although an 
excess can be seen in the ce ntral cell, which corresponds 
to the overdensity detected bv lFroebrich. Scholz fc Raftervl 
(|2007t ). Cz23, on the other hand, clearly detaches in the 
central cell (top-right panel) against a less-irregular field. 
The bottom panels show the surface densities built with 
field-star decontaminated photometry (Sect.[3]). 

4. Field-star decontamination 

Especially in the case of FSR 1716, the stellar surface- 
density in the direction of both objects clearly shows that 
field-star contamination should be taken into account. This 
assertion is further supported by the J x ( J — H) and 
J X {J — Ks) CMDs extracted from the R < 3' region of 
FSR 1716 (Fig. [3D, and the R < 4' region of Cz23 (Fig. [J). 
Features present in the central CMDs and in the respective 
comparison field (top and middle panels), show that field 
stars contribute in varying proportions to the CMDs, espe- 
cially for the bulge-projected FSR 1716. Nevertheless, when 
contrasted with the equal-area comparison field extractions, 
sequences typical of old stellar systems are suggested in 
both cases, a conspicuous giant branch in FSR 1716 and a 
giant clump in Cz23. 

Field-star decontamination is a very important, yet dif- 
ficult, step in the identification and characterisation of 
star clusters. Different approach es have been u sed to this 
purpose, among the m, those of iMercer et al.l (|2005D and 
ICarraro et al.l (|2006f) . The first is based essentially on spa- 
tial variations of the star-count density, but does not take 
into account CMD properties. In the latter, stars in a CMD 
extracted from an assumed cluster region are subtracted ac- 



- Observed (R=3') 




g Field Stars (same area) 




0.0 0.3 0.6 0.9 1.2 1.5 0.2 0.6 1.0 1.4 1.8 2.2 



(J-H) (J-K3) 

Fig. 3. 2MASS CMDs extracted from the i? < 3' region 
of FSR 1716. Top panels: observed photometry with the 
colours Jx{J -H) (left) and Jx{J - Ks) (right). Middle: 
equal-area (19.77' < R < 20') extraction of the comparison 
field, where the important disk and bulge contamination 
can be seen. Bottom panels: decontaminated CMDs that 
reveal a highly reddened, relatively populous MS and a gi- 
ant branch, typical of old clusters. 

cording to colour and magnitude similarity with the stars 
of an equal-area comparison field CMD. 

In the p resent case, we apply th e statistical algorithm 
described in iBonatto fc Bical (|2007b( ) to quantify the field- 
star contamination in the CMDs. The algorithm makes 
use of both approaches above, in the sense that relative 
star-count density together with colour/magnitude similar- 
ity between cluster and comparison field extractions are 
taken into account simultaneously. It measures the rela- 
tive number densities of probable field and cluster stars 
in cubic CMD cells whose axes correspond to the J mag- 
nitude and the (J — H) and (J — Kg) colours. These are 
the 2MASS colours that provide the maximum variance 
am ong CMD seq uences for OCs of different ages (e.g. 
Bo natto. Bica fc G irardi 2004). The algorithm: (i) divides 
the full range of magnitude and colours covered by the 
CMD into a 3D grid, (ii) calculates the expected num- 
ber density of field stars in each cell based on the num- 
ber of comparison field stars with similar magnitude and 
colours as those in the cell, and (Hi) subtracts the ex- 
pected number of field stars from each cell. The algorithm 
is responsive to local variations of field-star contamina- 




Fig. 4. Same as Fig. [3] for the region i? < 4' of Cz 23. 

Contamination, in this case, is less important. The equal- 
area comparison field extraction was taken from the region 
19.6' < R < 20'. A giant clump shows up especially in the 
decontaminated CMDs, denoting advanced age. 

tion (jBonatto fc Bicall20071j ). Typical cell dimensions are 
AJ = 1.0, and A(J - H) = A(J - Ks) = 0.25, which are 
wide enough to allow for sufficient star-count statistics in 
individual cells and narrow enough to preserve the morphol- 
ogy of the CMD evolutionary sequences. The comparison 
fields are located within R = 15' - 30' (FSR1716), and 
R = 10' — 30' (Cz 23). In both cases, the inner boundary of 
the comparison field lies beyond the probable tidal radius 
(Sect, in]), which minimises the probability of oversubtrac- 
tion of member stars. We emphasise that the equal- area 
field extractions shown in the middle panels of Figs. [3] and 
m serve only for comparisons among the panels. Actually, 
the decontamination process is carried out with the wide 
surrounding area as described above. Further details on the 
algorithm, including discuss ions on subtract i on effic iency 
and limitations, are given in 'Bonatto & Bica' f2007 bD. 

As extensively discussed in Bonatto & Bica ( 2007U) . 
differential reddening between cluster and field stars is re- 
ally critical for the decontamination algorithm. Large gradi- 
ents would require wide cell sizes or, in extreme cases, would 
preclude application of the algorithm altogether. Basically, 
it would be required, e.g. |A(J — H)\ cell size (0.25, in 
the present work) between cluster and comparison field for 
the differential reddening to affect the subtraction in a given 



Fig. 5. Same as Fig. [3]for the i? < 3' region of a 'test-field' 
located at 1° to the East of FSR 1716. The decontaminated 
CMDs (bottom panels) contain essentially statistical noise. 

cell. However, in both cases the CMDs extracted from the 
cluster region and comparison field (Figs. [3] and [3]) indicate 
that the differential reddening, if present, is not important. 

The decontaminated stellar surface-density distribu- 
tions (bottom panels of Fig. [2]) are an indicator of the algo- 
rithm efficiency. In both cases the central cells, which corre- 
spond to the location of the overdensities, present conspic- 
uous surface densities. Besides, the residual surface-density 
around the centre was reduced to a minimum level. We note 
that the decontamination is essentially based on the colour- 
magnitude distributions of the stars in different spatial re- 
gions. The fact that the decontaminated surface density 
ends up with a conspicuous excess only at the assumed clus- 
ter region (Fig. [2|) implies significant differences among the 
spatial regions, both in the colour-magnitude distributions 
and in the number of stars within a given colour-magnitude 
bin. This is expected of star clusters, which can be basi- 
cally characterised by a single-stellar population, projected 
against a Galactic stellar field. 

The decontaminated CMDs are shown in the bottom 
panels of Figs. [3] and [H As expected, essentially all of 
the bulge and disk contamination in FSR 1716 is removed, 
leaving stellar sequences typical of a reddened old OC, 
with a rather populous main sequence (MS), well-defined 
main sequence turnoff (MSTO), and the giant branch. 
Alternatively, the decontaminated morphology might re- 
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semble that of a poorly-populated GC. We explore this 
possibility in Sect. 15.11 Of the 588 stars present in the 
i? < 3' CMD, only 128 remain. A similar conclusion ap- 
plies to Cz23, in which the disk contamination has been 
subtracted, reveahng a poorly-populated, old OC; 66 of the 
228 stars remain in the CMD after decontamination. 

The approximately central direction of FSR 1716 and 
the poorly-populated nature of Cz 23 require additional sta- 
tistical analysis. For this purpose, we present in Table[2]the 
full statistics of the decontaminated sequences and field 
stars, discriminated by magnitude bins. Statistically rel- 
evant parameters that characterise the nature of a star 
cluster are: (i) Niu which, for a given magnitude bin, 
corresponds to the ratio of the decontaminated number 
of stars to the la Poisson fluctuation of the number of 
observed stars, (ii) ctfSi which is related to the proba- 
bility that the decontaminated stars result from the nor- 
mal star count fluctuation in the comparison fleld and, 
(iii) FSunifi which measures the star-count uniformity of 
the comparison fleld. Properties of Nia-, <tfSj and FSunif, 
measured in OCs and flel d fluctuations are discussed in 
iBica. Bonatto fc Camargd ( 



. Table [2] also provides in- 
tegrated values of the above parameters, which correspond 
to the full magnitude range spanned by the CMD of each 
OC. The spatial regions considered here are those sampled 
by the CMDs shown in the top panels of Figs. [2] and IH 

Star cluster CMDs should have integrated A^i,^ val- 
ues s ignificantly higher than 1 (jBica. Bonatto fc Camargol 
I2008D . a condition that is met by FSR 1716 {Ni^ = 5.0) and 
Cz23 (iViCT = 4.4). Besides, the number of decontaminated 
stars in each magnitude bin of FSR 1716 is higher (at the So- 
level or higher) than what could be expected from field-star 
fluctuations. As a further test of the statistical signiflcance 
of the above results we investigate star count properties of 
the fleld stars. First, the comparison field is divided into 
8 sectors around the cluster centre. Next, we compute the 
parameter dps ; which is the 1 a Poisson fluctuation around 
the mean of the star counts measured in the 8 sectors of 
the comparison field (corrected for the different areas of the 
sectors and cluster extraction) . In a spatially uniform com- 
parison fleld, (Tps is expected to be low. In this context, star 
clusters should have the probable number of member stars 
(A^ci) higher than ~ 3 crps, to minimise the probability that 
A^ci arises from fluctuations of a non-uniform comparison 
field. This condition is fully satisfied, in some cases reach- 
ing the level A'ci ~ 5 crps ■ The ratio decreases somewhat 
for Cz 23, probably because it has a low number of mem- 
ber stars. Finally, we also provide in Table[2]the parameter 
FSunif- For a given magnitude bin we first compute the 
average number of stars over all sectors (N) and the cor- 
responding 1(7 fluctuation ct^tv); thus, FSunii is deflned as 
FSunif = cr^N) / (N) . Non uniformities such as heavy differ- 
ential reddening should result in high values of F^unif • 

Since we usually work with comparison fields wider than 
the possible-cluster extractions, the correction for the dif- 
ferent spatial areas between field and cluster is expected 
to produce a fractional number of probable field stars 
(n^g'^g^p) in some cells. Before the cell-by-cell subtraction, 
the fractional numbers are rounded off to the nearest in- 
teger, but limited to the number of observed stars in each 
cell n^'^"^^^ = ^ H'n'YsLxp) - ^ots^ where NI represents the 
nearest integer. The global effect is quantified by means of 
the difference between the expected number of field stars 
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Fig. 6. Left panels: Observed and decontaminated CMDs 
extracted from R — 4' comparison fields, taken at 20' to 
each side of Cz 23. The corresponding RDPs (right panels) 
exhibit essentially field-fiuctuation features. 



in each cell {nYslexp) actual number of subtracted 

stars (nYs^sub)- Summed over all cells, this quantity pro- 
vides an estimate of the total subtraction efficiency of the 
process. 



fsub = 100 X 



cell 



cell 



cell 

'fs^exp 



(%). 



Ideally, the best results would be obtained for an efficiency 
fsub ~ 100%. With the assumed grid settings for the de- 
contamination of FSR 1716 and Cz23, the subtraction effi- 
ciencies turned out to be higher than 93%. 

Finally, it's worth noting that the qualitative and quan- 
titative expectations of the decontamination algorithm 
have been met by the output. On the one hand, the de- 
contaminated photometry reveals a very-high excess in the 
surface-density distributions with respect to the surround- 
ings, in both cases (Fig. [2]). In addition, CMDs extracted 
from the spatial regions where the excesses occur present 
statistical significant (Table [2]) cluster- like features (Figs. [3] 
andlH). 

4.1. Additional test for FSR1716: field at 1° to the East 

Evidence drawn from the previous sections indicates that 
FSR 1716 is an old star cluster. However, since it lies at 
a low latitude and is projected not far from the dense 
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Table 2. Statistics of the field-star decontamination discriminated by magnitude bins 
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Table Notes. The table provides, for each magnitude bin (AJ), the number of observed stars (A''obs) within the spatial region 
sampled in the CMDs shown in the top panels of Figs. [3] and |4j the respective number of probable member stars (A^d) 
computed by the decontamination algorithm, the Ni^ parameter, the la Poisson fluctuation (crps) around the mean, with 
respect to the star counts measured in the 8 sectors of the comparison fleld, and the fleld-star uniformity parameter. The 
statistical signiflcance of A^ci is reflected in its ratio with the la Poisson fluctuation of A'obs (A'ict) and with aps- The bottom 
line corresponds to the full magnitude range. 



bulge stellar field, we provide an additional test to fur- 
ther probe its nature. A 'test-field' of 20' of radius was 
extracted at 1° to the East of FSR 1716, with Galactic 
coordinates £ = 330.79° and b = —1.59°, thus projected 
somewhat closer to the centre than FSR 1716. The 'test- 
field' photometry was analysed in the same way as that 
of FSR 1716. Similarly to FSR 1716, we considered a cen- 
tral region (i? < 3') of the 'test-field', extracted the equal- 
area comparison field near the border (19.77' — 20'), and 
applied the decontamination algorithm. The results are 
shown in Fig. [5l Contrarily to the old-cluster decontam- 
inated CMD morphology of FSR 1716 (Fig. E]), the de- 
contaminated CMD of the 'test-field' (bottom panels of 
Fig. ^ essentially contains Poisson noise produced by the 
mutual subtraction of star samples that share similar dis- 
tributions of magnitude and colours. The same conclusion 
can be drawn from the featureless 'test-field' RDF (Sect. [SI 
Fig. [101). 

As discussed in iBica. Bonatto fc Camargol (j2008l) . 
Poisson fluctuations of the dense stellar field projected to- 
wards the bulge do not produce cluster-like CMDs and 
RDPs simultaneously. In particular, when field-fluctuation 
CMDs are field-star decontaminated, what results is essen- 
tially represented by the above FSR 1716 'test-field' exper- 
iment (Fig. [5]). In this sense, the striking differences exhib- 
ited by the CMD and RDP of FSR 1716, as compared to 
those of the 'test-field', can be taken as a robust evidence 
of the old-cluster nature of FSR 1716. 



4.2. Additional comparison fields for Cz23 

The CMD of Cz23, especially the decontaminated one 
(Fig.[4]), suggests a poorly-populated old OC. Thus, we test 
whether the decontamination procedure, applied to ran- 
dom fields around Cz23, may produce similar cluster- like 
sequences and stellar radial profiles. We consider R — 4' 
comparison fields, taken at 20' to each side of Cz 23. These 
fields have the same projected area as that of Cz 23 shown 
in Fig. [H and the centre to centre offsets correspond to « 4 
times the cluster radius of Cz23 (Table [3]). Field-star de- 



contamination was applied to these CMDs using the same 
offset field as that employed for Cz23 (Sect. [J]), but re- 
stricted to i? = 25' — 30' to avoid self subtraction. 

The observed and decontaminated comparison field 
CMDs are shown in Fig. [Hi In stark contrast with the de- 
contaminated CMD of Cz23 (Fig. [3]), what remains in the 
R — A' comparison field CMDs is a randomly scattered 
small number of stars, similar to the CMD of the FSR 1716 
'test-field' (Sect. ITT|) . To complete this analysis we also in- 
clude in Fig. [Sj the RDPs of the comparison fields. Both 
CMDs and RDPs of the comparison fields are typical of 
Poisson fluctuations of stellar fields projecte d towards the 
Galactic anti-centre ('Bon atto fc Bicall2008bD . 

As expected of field fluctuations, the CMDs and RDPs 
of the comparison fields (Fig. [S]) are featureless, quite con- 
trasting with the cluster CMD (Fig. [3]) and RDP (Fig. [11]) 
of Cz23. 

5. Age, reddening, and distance 

Having established that FSR 1716 and Cz 23 present old- 
age features, we proceed to derive their fun damental pa- 
rame ters. We work with Padova isochrones (jCirardi et al.l 
|2002[ ) of Solar and sub-Solar metallicity, computed with 
the 2MASS J, H, and filterfl These isochrones 
are very similar to the Johnson-Kron-Cousins ones (e.g. 
JBcssel & Brett 1988), with differences o f at most 0.01 in 
(J - H) (Bonatto, Bica & Gii^M\2QQi)- We adopt Rq = 
7.2±0.3kpc (.Bica et al.n2006i) as the Sun's distance to the 
Galactic centre to compute the OCs Galactocentric dis- 
tances. The value Rq = 7.2±0.3kpc was derived by means 
of the GC spatial distribution. Other re cent studies gave 
simil ar results, e.g. Rq = 7.2 ± 0.9 kpc (lEisenhauer et al.l 
12003? ), Rq = 7.62 ± 0. 32 kpc (lETsenhauer et al.ll2005[ ) and 
Rq = 7.52 ± 0.10 kpc ([Nishivama et al.ll2006[ ). with differ- 
ent approaches. We make fits by eye, taking into account 
that, because of the presence of binaries, the isochrone 
should be shifted to the left of the MS fiducial fine (e.g. 
iBonatto. Bica fc Santos Jrll2005l and references therein). 



http:/ /stev. oapd. inaf. it/cgi-bin/cm(l\ 
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Fig. 7. Top panels: the field-decontaminated {R < 3') 
CMD morphology of FSR 1716 (empty chxles) is com- 
pared to that of the ~ 7Gyr OC NGC 188 (gray circles). 
Jx (J - H) (left panels) and Jx (J - K^) (right) CMDs are 
shown. Bottom: same as above for the i? < 4' CMD of Cz 23 
(empty circles) and M67 (gray circles). The adopted fit 
(Sect. [5XT|) for FSR 1716 corresponds to the 7Gyr, Solar- 
metallicity Padova isochrone, while for Cz23, we adopted 
the 5Gyr, Solar-metallicity isochrone (Sect. I5.3.2p . The 
shaded polygon overplotted in the left panels shows the 
colour-magnitude filter. 

In principle, the decontaminated CMD morphologies 
(bottom panels of Figs. [3] and 2]) - as well as the additional 
tests discussed in Sect. 14.11 - should provide enough con- 
straints to derive reliable cluster fundamental parameters. 
However, in view of the importance of identifying new old 
star clusters, we explore a range of possibili ties i n terrn s 
of age and metalhcity. As discussed by, e.g. iFriell (|1995f) . 
OC metallicities in general range from Solar {[Fe/H] — 0, 
OT Z = 0.019) to sub-Solar {[Fe/H] = -0.5, Z = 0.006) 
values. 

5. 1 . CMD morphology of FSR 1716 

First, we compare in Fig. [7] (top panels) the decontaminated 
CMD of FSR 1716 to that of the - 7Gyr OC NGC 188 
(|Bonatto. Bica fc Santos Jr]|2005l and references therein). 
CMD morphologies in both colours show a good agreement 
along the MS, MSTO, giant branch, and the MS width, 
which suggests that both clusters are of similar age. In 
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Fig. 8. Age and metallicity interpretation of FSR 1716. 
Top-left panel: isochrones of different ages and fixed Solar 
metallicity. Top-right: same ages but with the sub-Solar 
metalhcity [Fe/H] = -0.5 {Z = 0.006). Fits have in 
common a similar description of the MS. Middle pan- 
els: template fit assuming a GC-like RGB and sub-giant 
branch morphology for FSR 1716, with the GC NGC 6397 
{[Fe/H] = -1.95, Z = 0.0002 - left panel), and M4 
{[Fe/H] = -1.21, Z = 0.0012 - right). Padova isochrones of 
12 Gyr and same metallicity as the GCs are shown. Bottom 
panels: same as above but assuming a brighter MSTO for 
FSR 1716. The colour range in each panel is optimised to 
shown differences among the fits (top) and enhance the rel- 
evant features (bottom). 



particular, both CMDs present a similar smoothly-curved 
MSTO and a scarcely populated red clump. 

Different ages and metallicities are tested in Fig. [8l 
Solar-metallicity solutions with isochrones of ages 5, 7, and 
9 Gyr are shown in the top-left panel. Same ages but with 
the low OC metallicity range {[Fe/H] ^ —0.5, which cor- 
responds to Z « 0.006), are considered in the top-right 
panels. As fit constraint, the isochrones have in common a 
similar description of the MS. In both panels, the youngest 
isochrone produces a relatively poorer fit, especially to the 
giant branch. On the other hand, the older isochrones - of 
both metallicity ranges - produce qualitatively similar fits. 
This not only confirms the old age of FSR 1716, but, to- 
gether with a few blue stars at J ~ 14.3 and (J — H) « 0.3 
(Fig. [7]), which might suggest a blue horizontal branch 
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Fig. 9. Age and metallicity interpretation of Cz23. Top 
panels: isochrones of different ages and fixed Solar metal- 
licity. Bottom: same ages but with the sub-Solar metallicity 
[Fe/H] = -0.5 {Z = 0.006). Fits are required to provide a 
similar description of the MS. 



(HB), raise the possibility of a poorly-populated, metal- 
poor GC. 

To test the GC hypothesis we compare the decontam- 
inated CMD of FSR 1716 to those of the blue HB GCs 
NGC6397 and M4 (NGC6121). Both GCs are among the 
nearest ones (rf© « 2.3 kpc), which provides MS depth for 
2MASS photometry, and are metal-poor, with [Fe/H] ^ 

-1.95 (NGC6397) and [Fe/H] 1.2 (M4). We also use 

Pado va isochrones of 12 G yr to characterise the GC ages 
(e.g. lOrtolani et a"nil995D . The CMDs of NGC6397 (ex- 
tracted within i?ext < 3') and M4 (i?cxt < 2') were pro- 
duced in the same way as those of FSR 1716 (Sect. [3]). 
For a quantitative derivation of fundamental parameters by 
means of a comparison with GCs, the isochrones were first 
set to the GCs according to the fundamental data taken 
from Harris (1996, and the 2003 updat^H). Subsequently, 
we searched for the best overall match between the GC -f 
isochrone and FSR 1716 sequences. The results are shown 
in the middle panels of Fig. [8l As a further test, we also 
searched for solutions with the same isochrones and tem- 
plate GCs as above, but assuming a brighter MSTO for 
FSR 1716, located at J « 13.5. These tentative fits are sown 



in the bottom panels of Fig. [ffl. Compared to NGC 6397 
(bottom- left panel), the general fit is acceptable, although 
the RGB of FSR 1716 is somewhat redder, which might 
suggest a higher metallicity. However, the fit with the 
more metal-rich GC M 4 (bottom-right) is obviously poorer. 
Besides the red giants which are too red as compared to 
M4, the sub-giant branch morphologies become critically 
different. In the case of the M4-like solution, FSR 1716 
would be as close (d© « 0.6 kpc) as an old OC. In the case 
of the NGC 6397 solution, FSR 1716 would be the nearest 
(^0 ~ 0.5 kpc) GC so far discovered. 

The above comparisons show that a very low metal- 
licity, like that of NGC 6397, is not compatible with the 
giant branch of FSR 1716, since it is steeper and bluer in 
NGC 6397. M4, on the other hand, provides a better match 
for the giant branch. What makes this possibility less likely 
is that the few extreme blue stars in the CMD of FSR 1716 
do not have counterpart in either NGC 6397 or M 4, whose 
blue HBs are redder and well distributed in colour. The 
CMD shape for 13 ^ J ^ 15 in FSR 1716 is more com- 
patible with the MSTO and MS of an old OC than the 
matched subgiant branches of the two GCs. However, the 
evidence of a Palomar-like (i.e. l ow-mass and loose struc- 
ture) GC (jBonatto &: Bicair2008a[ ) , although marginal, can- 
not be ruled out with the existing data. 

We call attention to the fact that it is not an easy task 
to establish the nature of a low-mass GC or a low-mass 
old OC based on a few candidate HB or giant clump stars. 
For instance, NGC 188 (and other old O Cs) presents only 
a few UV-bright stars, as shown by, e.g. Landsman et al.l 
()l998l ). In fact, both low-mass GCs and low-mass old 
OCs can hardly have any HBs or clump giants, includ- 
ing NGC 188 with only 5 clum p giants (iDinescu et al. I 
1995t and lLandsman et al.lll998( ). As examples of low-mass 
GCs with a small number o f UV-bright stars we quote 
P a.lomarl3 (ICote et al. 1120021) with only a few, and AM- 
4 (jlnman fc CarnevI 1987 ) with none. 

Given the above circumstances, deeper observations 
would be necessary for a more conclusive definition on the 
nature of FSR 1716 as a GC (as originally suggested by 
iFroebrich. Scholz fc RaftCT^[200l or a very old OC, with 
a metallicity that can be sub Solar. 



5.2. CMD morphology of Cz23 

The CMD features of Cz23 (bottom panels of Fig. [T]) 
indicate a younger age than FSR 1716. Indeed, its field- 
decontaminated CMD morphology is similar to that of the 
old OC M67. The age of M 67, as derived from ne ar- infrared 
photometry, is ~ 3.2 Gyr (jBonatto fc Bicall2003l . and refer- 
ences therein). More recent determinations based on spec- 
troscopy of individual stars put the age at ^ 4 Gyr (e.g. 
iGiampapa eralll2006[ ). In addition, the MSTO of Cz23 
seems to be at a slightly fainter magnitude range, which 
indicates an age older than that of M 67. 

Thus, to probe the age and metallicity of Cz23, we 
show in Fig. [9] fits with the 3, 5, and 6 Gyr isochrones 
with a fixed Solar metallicity (top panels), and the sub- 



http:/ /physun. physics. mcmaster.ca/GlobuIar. html 



Although at approximately the same distance from the Sun 
and absorbed by a similar reddening value, the 2MASS cutoffs of 
M 4 and NGC 6397 are somewhat different. This occurs basically 
because M4 is projected against the bulge {£ « 351°), while 
NGC 6397 lies at the outskirts {£ « 338°). 
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Solar metallicity [Fe/H] ^ —0.5 (bottom). The fits are 
required to provide a similar description of the MS. It is 
clear that, under this condition, the 3Gyr isochrone fails 
to reproduce the giant clump, in both metallicity ranges. 
The 5 and 6Gyr isochrones, on the other hand, provide ac- 
ceptable MS and clump fits, which confirms the old age of 
Cz 23. The 6 Gyr solutions are obtained with an essentially 
zero reddening applied to the isochrones, which implies that 
older isochrones would require negative correction. This ar- 
gument constrains the age of Cz23 to the range « 4 to 
« 6 Gyr. 

5.3. The adopted age/metallicity solutions 

5.3.1. FSR 1716 

Most of the arguments drawn above favour the old OC 
interpretation for the nature of FSR 1716. Since metal-poor 
OCs are p referentially distributed in the outer Galaxy (e.g. 
Iaie][i99i), and FSR 1716 is located inside the Solar circle, 
its probable metallicity should be around the Solar value. 

Based on the above, we take the 7 Gyr, Solar metal- 
licity isochrone as representative of the stellar sequences 
of FSR 1716, with an uncertainty of about ±lGyr. The 
corresponding best-fit was obtained with the reddening 
E{J — H) — 0.63 ± 0.02, and observed distance modulus 
(m — M)j = 11.1 ± 0.2. This solution is shown in the top 
panels of Fig. [T) Taking into account fit uncertainties we de- 
rive the age 7 ± 1 Gyr, reddening E{B -V)= 2.0 ± 0.1 or, 
equivalently (Sect. [3]) Ay — 6.3 ±0.2, the absolute distance 
modulus {m — M)o = 9.36 ±0.21, and the distance from 
the Sun dg = 0.8 ± 0.1 kpc. Thus, for Rq = 7.2 kpc, the 
Galactocentric distance of FSR1716 is Rqc = 6.6±0.2kpc, 
which puts it w 0.6 kpc inside the Solar circle. 

Alternatively, we consider as well the less-likely pos- 
sibility of a GC. In this case, the best-fit of the 12 Gyr 
isochrone (based on the similarity with the M 4 CMD mor- 
phology - Fig. [5]) is obtained with (m - M)j = 13.5 ± 0.3, 
E{J -H)^ 0.63 ± 0.04, E{B -V)= 2.0 ± 0.1 and Ay = 
6.3 ± 0.2, (m - M)o = 11.76 ± 0.32, = 2.3 ± 0.3 kpc, 
and i?GC = 5.4 ± 0.2 kpc, thus ~ 1.8 kpc inside the Solar 
circle. 

5.3.2. Cz23 

With similar arguments as those used to estimate the age 
of FSR 1716, we take the 5 Gyr (±lGyr) isochrone to 
represent the age of Cz 23. Since both metallicity ranges 

{[Fe/H] 0.5 and 0.0) produce similar CMD fits (Fig.©, 

for simplicity, we also work with the Solar metallicity 
isochrone. 

Thus, fundamental parameters of Cz 23 are E{J — H) = 
0.00 ± 0.01, E{B -V) = 0.00 ± 0.03 and Ay = 0.0 ± 0.1, 
(m - M)j = (m - M)o = 12.0 ± 0.1, = 2.5 ± 0.2 kpc, 
and i?GC = 9.7 ± 0.2 kpc; Cz23 hes « 2.5 kpc outside the 
Solar circle. This solution is shown in the bottom panels of 
Fig. El 



6. Cluster structure 

Structural parameters are derived by means of the RDPs, 
defined as the projected radial distribution of the number 
density of stars around the cluster centre. 



Star clusters usually have RDPs that follow a well- 
defined analytical profile. The most often used are the single 
mass, modified isot hermal spher e ((King 1966), the modified 
isoth ermal sphere (IWilsonlll975f). a nd the power law with a 
core (jElson. Fall fc FreemanI 1987t ). Each function is char- 
acterised by different parameters that are related to cluster 
structure. However, because the error bars in the present 
RDPs are significant (Fig. fTU)) . we use the analytical profile 
a{R) = (Tbg + o'o/(l + {R/RcY), where (Jhg is the residual 
background density, ctq is the central density of stars, and 
Rc is the co re radius. This function is similar to that by 
iKingj (|1962| ) to describe the surface brightness profiles in 
the central parts of globular clusters. 

The RDPs are built with colour-magnitude filtered pho- 
tometry to minimise noise. The most probable cluster se- 
quences are isolated by means of colour-magnitude fil- 
ters, which are used to exclude stars with colours dif- 
ferent from those of the assumed cluster sequence. They 
are wide enough to include cluster MS and evolved star 
colour distributions, as well as the la photometric un- 
certainties. Colour-magnitude filter widths should also ac- 
count for formation or dynamical evolution-related ef- 
fects, such as enhanced fractions of binaries (and other 
multiple systems) towards the central parts of clus- 
ter s, since suc h system s tend to widen the MS (e.g. 
Hurlev fc ToutI 11998: Kerber et al.] 120021: iBonatto fc Bical 

120051) . 



I2007bl : IBonatto. Bica fc Santos Jr.l l2005f) . The colour- 
magnitude filters for the present OCs are shown in the left 
panels of Fig. [71 Residual field stars with colours similar 
to those of the cluster are expected to remain inside the 
colour-magnitude filter region. They affect the intrinsic stel- 
lar RDP in a way that depends on the relative densities of 
field and cluster stars. The contribution of the residual con- 
tamination to the observed RDP is statistically subtracted 
by means of the field. In practical terms, the use of colour- 
magnitude filters in cluster sequences enhances the contrast 
of the RDP with respect to the background, espec ially for 
clusters in dense fields (e.g. IBonatto fc Bicall200T ^. 

Oversampling near the centre and undersampling at 
large radii are avoided by using rings of increasing width 
with distance from the centre. A typical set of ring widths 
is Ai? = 0.5, 1, 2, and 5', respectively for 0' < i? < 1', 
1' < i? < 4', 4' < i? < 10', and 10' < i? < 30'. The number 
and width of the rings can be set to produce RDPs with ad- 
equate spatial resolution and as low as possible la Poisson 
errors. The residual background level of each RDP corre- 
sponds to the average number of colour-magnitude filtered 
stars measured in the field. The R coordinate (and uncer- 
tainty) of each ring corresponds to the average position and 
standard deviation of the stars inside the ring. 

The RDPs of FSR 1716 and Cz 23 are given in Fig. [lOl 
Because of the uncertainties associated to the age (and 
metallicity) derivation (Sect. [5]), which propagate to the ab- 
solute value of the structural parameters, RDPs in Fig. [TOl 
are shown in angular scale. Besides the RDPs built with 
the colour-magnitude filters, we show, for illustrative pur- 
poses, those produced with the observed (raw) photome- 
try. Especially for FSR 1716, minimisation of the number 
of non-cluster stars by the colour-magnitude filter resulted 
in a RDP with a higher contrast with respect to the back- 
ground. Fits of the King-like profile were performed with 
a non-linear least-squares fit routine that uses errors as 
weights. To minimise degrees of freedom, ctq and Rc were 
derived from the RDP fit, while abg is measured in the 
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Table 3. Cluster structural parameters 



Cluster 


1' 






<5c 




Rc 


^^RDP 




-Rrdp 




(pc) 


(stars arcmin^^) 


(stars arcmin^^) 






(') 


(') 


(pc) 


(pc) 


(1) 


(2) 


(3) 


(4) 


(5) 




(6) 


(7) 


(8) 


(9) 


FSRiriet 


0.216 


10.5 ± 1.0 


16.2 ±4.8 


2.6 ±0.5 


0.8 


8 ±0.23 


6.0 ±0.6 


0.19 ±0.05 


1.3 ± 0.2 


FSR1716t 


0.653 


10.5 ± 1.0 


16.2 ±4.8 


2.6 ±0.5 


0.8 


8 ±0.23 


6.0 ±0.6 


0.57 ±0.15 


3.9 ±0.4 


Cz23 


0.728 


2.8 ±0.1 


16.5 ±4.7 


7.0 ± 1.7 


0.49 ± 0.04 


4.9 ±0.7 


0.36 ±0.08 


3.6 ±0.4 



Table Notes. Col. 2: arcmin to parsec scale. King profile is expressed as o{R) — atg+ao/il + iR/Rof). To minimise degrees of 
freedom in RDP fits, (Ji,g was kept fixed (measured in the respective comparison fields) while ctq and Rc were allowed to vary. 
Col. 5: cluster/background density contrast {Sc = 1 + ao/utg), measured in CM-filtered RDPs. Cols. 6-9, core and cluster radii 
in angular and absolute units. Parameters of FSR1716 are derived for the OC (f) or GC (J) interpretation (Sect. [5]). 



field. These values are given in Table [31 and the best-fit so- 
lutions are superimposed on the colour-magnitude filtered 
RDPs (Fig. (TU]) . Table [3] also provides structural parame- 
ters in absolute units, computed with the cluster distances 
derived in Sect. 15.31 Because of the 2MASS photometric 
limit, which in most cases corresponds to a cutoff for stars 
brighter than J « 16, cto should be taken as a lower limit 
to the actual central number density. The adopted King- 
like function describes well the RDPs along the full radius 
range, within uncertainties. 

To complete the structural description of the clusters 
we also estimate the cluster radius (i?RDp) and uncer- 
tainty by visually comparing the RDP level (and fluc- 
tuations) with the background. The cluster radius cor- 
responds to the distance from the cluster centre where 
RDP and background are statistically indistinguishable 
(e.g. iBonatto fc Bical 1200 5'. and references therein). For 
practical purposes, most of the cluster stars are contained 
within i?RDP- Note that i?RDP should not be mistaken 
for the tidal radius. Tidal radii are derived from King 
fits to RDPs, which depend on wide fields and adequate 
Poisson errors. For instance, in populous and relatively 
high Galactic latitude OCs such as M67, NGC 188, and 
NGC2477, cluster ra dii are a factor ^ 0.5 - 0.7 of the 
respective tidal radii (jBonatto fc Bical |2005[ ). The cluster 
radii of the present objects are given in cols. 6-7 (angular 
scale) and 8-9 (absolute scale) of Table [3l If cluster and 
tidal radii of FSR 1716 and Cz 23 are similarly related as 
in the bright OCs, the lower-limit of the radial range taken 
as comparison field (Sect. [4]) is located beyond the proba- 
ble tidal radius. This, in turn, minimises the probability of 
cluster members at large radii, e.g. in the cluster halo, to be 
considered as field stars by the decontamination algorithm. 

Compared to the distribution of core radius derived for a 
sample of relatively near by OCs bv lPiskunov et all (|2007l) . 
FSR 1716 (especially the OC interpretation) and Cz 23 oc- 
cupy the small- i?c tail. Assuming that the tidal radius is 
~ 2xi?RDP, Cz23 and FSR 1716 (GC interpretation) would 
be around the median value of the distribution, while the 
OC interpretation of FSR 1716 again would be at the small- 
tidal radius tail. 

Table|3](col. 5) provides the density contrast parameter 
6c = I +o'o/o'f)g. Since Sc is measured in colour-magnitude- 
filtered RDPs, it may not correspond to the visual con- 
trast produced by observed stellar distributions in images. 
FSR 1716 presents a low contrast in the 2MASS Kg image 
(Fig.[T]) but, because most of the non-cluster stars have been 
excluded by the colour-magnitude filter, the correspond- 
ing RDP has a relatively high density contrast, Sc ~ 2.6. 
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Fig. 10. Stellar RDPs (filled circles) of FSR 1716 (left pan- 
els) and Cz23 (right). RDPs built with the observed (raw) 
photometry are shown in the top panels, while those built 
with colour-magnitude filtered photometry are in the bot- 
tom panels. Solid line: best-fit King-like profile. Horizontal 
shaded polygon: offset field stellar background level. Shaded 
regions: la King fit uncertainty. The core and cluster radii 
are indicated. Inset of panel (a): RDP of the test-field of 
FSR 1716. Angular scale is used. 



Obviously, the high-contrast OC Cz23 (Figs. [T] and [TU)) is 
reflected in the high value of the density-contrast parameter 
Sc ~ 7.0. Interestingly, FSR 1716 and Cz23 are projected 
against almost opposite directions (Table [T|) . Accordingly, 
the background (including foreground) stellar contribution 
turns out to be « 11 stars arcmin"^ in the direction of 
FSR 1716, and « 2.8 stars arcmin"^ towards Cz23. 
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Finally, the RDP of the 'test-field' of FSR 1716 is shown 
in the inset of panel (a) of Fig. (TUl As expected from the 
CMD analysis (Sect. 14. ip . the RDP is uniform from the 
centre to the borders of the field, characterised essentially 
by Poisson fiuctuations. The featureless 'test-field' RDP is 
low and flat, as compared the cluster RDP of FSR 1716 
(panels a and c). 

7. Mass function and cluster mass estimate 

Cluster mass functions {4>{m) = ^) are built follow- 
ing the methods presented in iBonatto fc Bical (|2005( ) 
(and references therein). We build them with colour- 
magnitude filtered photometry, the 3 2MASS bands 
separately, and the mass-luminosity relations obtained 
from the Padova isochrones and distances from the 
Sun adopt ed for FSR 1716 (Sect. IKXTj) and Cz23 
(S ect. 15. 3. 21). Further details o i i MF construction are given 
in lBica. Bonatto fc Blumberd (|2006l ). The effective magni- 
tude range over which MFs are computed is that where clus- 
ters present an excess of stars over the field. In both cases 
it begins below the MSTO and ends at a faint magnitude 
limit brighter than that stipulated by the 2MASS complete- 
ness limit (Sect. [3]). The effective MS stellar mass ranges are 
0.83 < m{MQ) < 1.1 (FSR 1716), and 0.94 < m{MQ) < 1.2 
(Cz23). 

The MFs computed for the whole cluster region {R < 
Rkdp) are shown in Fig.[TlJ Both cases suggest important 
depletion of the low-mass content, suggesting advanced dy- 
namical evolution, especially in the case of FSR 1716. The 
drop in the MF of FSR 1716 begins at m ^ 0.98 Mq, which 
corresponds to J w 14.5, more than 1 mag brighter than 
the 2MASS 99.9% completeness limit (Sect.|31) at the cluster 
position. This suggests that the drop may be real, although 
we cannot rule out crowding to account for part of the im- 
portant MF drop, since the 2MASS completeness limits re- 
fer to the field. Bearing in mind this caveat, we fit the MF of 
FSR 1716 with a two-segment function, (j){m) oc m~'^^'^^\ 
with X — —9.4 ± 0.6 in the mass range m < 0.98 Mq, and 
X = 2.3 ± 0.4 for 0.98 < to(Mo) < 1.1. The high-mass 
range slope is, with in the uncertai nty, somewhat steeper 
than the x — 1-35 of ISalpet"e^ (|l955l ) Initial Mass Function 
(IMF). However, the low-mass range slope is much flatter 
than Salpeter's IMF. The drop in the MF of FSR 1716 (at 
m « 1 Mq) a grees with one of the breaks present in the uni- 
versal IMF of lKroupal(|200H ). which assumes increasing flat- 
tening towards low-mass stars. This IMF is described by the 
slopes X = 0.3 ±0.5 for the range 0.08 < m(Afo) < 0.5 and 
X = 1.3 ± 0.3 for 0.5 < m(M0) < 1.0. The low-mass range 
of the MF of FSR 1716 is flatter than Kroupa's IMF, which 
again points to advanced dynamical evolution, crowding or, 
more probably, a combination of both. 

Although flat, the MF of Cz 23 is more monotonic than 
that of FSR 1716, and can be described by a single power- 
law, with the slope x = —3.1 ± 0.5, also flatter than 
Salpeter's or Kroupa's IMF. 

The number of observed MS and evolved stars in 
FSR 1716 is derived by counting the number of stars (in the 
background-subtracted colour-magnitude filtered photom- 
etry) that show up in the corresponding magnitude ranges, 
13.4 < J < 15.4 for the MS and J < 13.4 for the evolved 
stars. There are ums = 99 ± 5 and rievoi = 53 ± 7, MS 
and evolved stars, respectively; the corresponding mass val- 
ues are toms = 97 ± 5 Mq and nievoi = 58 ± 8 Mq . The 
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Fig. 11. Mass function computed for the whole region 
of FSR 1716 (top panel) and Cz23 (bottom). The MF 
of FSR 1716 can be represented by the broken power-law 
(j>{m) oc m~(i+x), with x = -9-4 ± 0.6 for m < 0.98 Mg, 
and X = 2.3 ± 0.4 for 0.98 < m{MQ) < 1.1. The MF drop 
begins at to « 0.98 Mq, which corresponds to J w 14.5. 
The MF of Cz23 can be described by a power-law with 
slope X — —3.1 ± 0.5. 



evolved star mass corresponds to rieuoi multiplied by the 
stellar mass at the MSTO, toto = 1-1 Mq. To estimate the 
total stellar mass we extrapolate the observed MF down to 
the H-burning mass limit (0.08 Mq). Because of the very 
fiat slope at the low-mass range, the extrapolated values 
are similar to the observed ones, Uextr = 160 ± 8 and 

niextr = 160±9Mq. 

The number of MS and evolved stars in Cz 23 is ums = 
54 ± 4 and rievoi = 5 ± 3, and the respective masses are 
nT-MS — 60±4Mq and nievoi = 7±4Mq; the extrapolated 
values are Uextr — 120 ± 8 and ruextr = 115 ± 7 Mq. 

Presently, both FSR 1716 and Cz 23 appear to be low- 
mass OCs, even when stars less massive than the observed 
range are taken into account. Most of the original mass 
must have been lost to the field because of dynamical effects 
along the several Gyrs since their formation (Sect. [8]). 

8. Discussion 

The range of acceptable isochrone solutions (Sect. [5]), 
together with the statistical tests (Sect. SJ applied to 
FSR 1716 and Cz 23, i.e. (i) the decontamination algorithm, 
(ii) the integrated and per magnitude iVio- parameter, and 
(Hi) the ratio of Nd to crpsi produce results consistent with 
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Fig. 12. Relations involving OC structural and fundamen- 
tal parameters. Circles: nearby OCs. Dotted circles: massive 
OCs (> IOOOM0). Two sets of parameters are shown for 
FSR 1716, according to the old OC and GC solutions. 



both objects being old OCs. Examination of photomet- 
ric and structural properties of the offset 'test-fields' for 
FSR 1716 and Cz23, also leads to the same conclusion. 

With the analyses of the preceding sections, we provide 
fundamental and structural parameters, most of which for 
the first time for FSR 1716 and Cz23. We now use these 
parameters to put both clusters into perspective, by com- 
paring some of their properties with those of a set of well- 
studied OCs. 

We compare the structural parameters, computed with 
the adopted isochrone solutions for FSR 1716 fSect. [573TT|) 
and Cz23 (Sect. I5.3.2p . with those of a reference sam- 
ple of nearby OCs with ages in th e range 70 — 7 OOP My r 
and masses within 400 - 5 300 (jBonatto fc Bicall2005b . 
To the original r eference sample were added the young 
OCs NGC 6611 (iBonatto Santos Jr. fc Bical I2006D and 
NGC4755 (jBonatto et al.1 l2006b[ ). Clusters are distin- 
guished ac cording to total mass (lower or higher than 
1 OOP M^). iBonatto fc Bical (|2005t) discuss parameter cor- 
relations in the reference sample. For completeness, both 
sets of parameters for FSR 1716 (Table [3]), corresponding 
to the old OC and GC possibilities, are considered sepa- 
rately. 

Core and cluster radii of the reference sample relate as 
^RDP = (8.9 ± P. 3) X ^(l o^o i) (panel a), which suggests a 
similar scaling for both kinds of radii, at least for the sam- 
pled ranges of P. 25 ^ Rc{pc) ^ 1.5 and 2 <, iiRDp(pc) ^ 15. 



Within uncertainties, FSR 1716 (both interpretations) and 
Cz 23 fit tightly into that relation. 

The reference OCs appear to follow the trend of increas- 
ing cluster radii with Galactocentric dis tance (panel b), a 
dependence previously suggested by, e.g. lLvng"ir(|l982l ). To 
explai n the increase of GC radii with Galacto centric dis- 
tance, Ivan den Be rgh. M orbev fc Pazdeil (|1991h suggested 
that part of the relation may be primordial, in the sense 
that the higher molecular gas density in central Galactic 
regions may have produced clusters with small radii. After 
formation, mass loss associated to stellar and dynamical 
evolution (such as mass segregation and evaporation), to- 
gether with tidal interactions with the Galactic poten- 
tial and giant molecular clouds, also contribute to the 
depletion of star clusters, especially the low-mass and 
centrally located ones (Sect. [T|). The cluster radius of 
FSR 1716 (GC interpretation) is consistent with the as- 
sumed relation. On the other hand, Cz23, and the OC 
interpretation of FSR 1716, appear exceedingly small for 
their Galactocentric distance, which is consistent with 
the several Gyrs of depletion. A similar dependence with 
Galactocentric distance holds as well for i?c, because of the 
relation with i?RDP implied by panel (a). 

In panels (c) and (d) of Fig. [T2] we compare the cluster 
and core radii of FSR 1716 and Cz 23 with those of the refer- 
ence sample in terms of age. The locus occupied by Cz 23 in 
both panels is consistent with the corresponding radii mea- 
sured in low-mass, old OCs. A similar conclusion applies to 
the GC interpretation of FSR 1716 (which again would be 
consis tent with t he structural parameters of Palomar-like 
GCs - IBonatto fc Bica 2PP8all . The significantly small radii 
of the OC interpretation of FSR 1716, especially 7?rdp, are 
consistent with those of OCs that h ave suffered severe de- 
pletion effects for long periods (e.g. IBonatto "fc'B ica 2PP5| 
IBonatto fc Bical[2PP7al and references therein). 

Another evidence in favour of dynamical effects af- 
fecting cluster size comes from the positions of FSR 1716 
and Cz23 in the Galactic plane fFig. [T51). Milky Way's 
spir al arm structure is based on iMomanv et al.l (j2006t ) 
and iDrimmel fc Spergell (|20(31[). der i ved fr om HII regions 
and molecular clouds (e.g. iRusseili 120031 ). The Galactic 
bar is shown with a 14° orien tation and 3 kpc in length 
(|Freudenreichl |l998: Vallcc'2PP5). For comparison purposes, 
we also include WEBDA clusters younger and older than 
1 Gyr. As discussed in Sect. [I] old OCs distribute preferen- 
tially outside the Solar circle. 

Interestingly, both clusters occur close to spiral arms, 
Cz 23 on Perseus, and FSR 1716 on Carina-Sagittarius (OC 
interpretation) or Crux-Scutum (GC). Since they are lo- 
cated close to the plane (TablejT]), they may have interacted 
with the arms, especially by means of encounters with gi- 
ant molecular cloudfl As discussed above, part of the small 
sizes of FSR 1716 and Cz23 ca n be accounted for by coUi 
sio ns with such structures (e.g. Wielen 19711: Wielen 1991 



Gi ejes et al.ll2PP6l : iGieles. Athanassoula fc Portegies Zwartj 
2PP^ . A similar effect was recently observed to occur with 
some small O Cs located close to th e Local and Orion- 
Cygnus Arms (jBonatto fc Bicall2PP8bD . 

Finally, both objects have MFs that appear to be very 
eroded and flat which, given the discussion on structural 
parameter above, reflects the consequences of several Gyrs 



Molecular clouds more massive than ^ lOr M(T) are found in 
the Solar neighbourhood (e.g. ISolomon et al.lll987l ). 
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Fig. 13. Galactic position of FSR 1716 and Cz 23 compared 
to the WEBDA OCs with ages younger (shaded circles) and 
older (filled circles) than 1 Gyr. The schematic projection 
of the Galaxy is seen from the North pole, with 7.2 kpc as 
the Suns distance to the Galactic centre. Main structures 
are identified. 



of the relentless dynamical effects and mass loss due to stel- 
lar evolution. Besides, both objects are low-mass clusters, 
with less than ~ 200 presently stored in stars. They 
must have been formed as more massive open clusters to 
have survived for so long in the Galaxy. 

9. Concluding remarks 

In this paper we study the nature of two stellar over- 
densiti es included in the catalogue o f can didate star clus- 
ters of iFroebrich. Scholz fc Raftervl (|2007D . FSR 1716 and 
FSR 834. Th e former was suggested to be a globu lar cluster 
candidate by IFroebrich. Scholz fc Raftervl ()2007( ). while the 
latter has the OC Cz23 as optical counterpart. The analy- 
ses are based on field-star decontaminated 2MASS CMDs 
and stellar radial density profiles, with algorithms previ- 
ously constructed by our group. Fundamental and struc- 
tural parameters of the clusters are derived. 

We present consistent evidence (e.g. CMD morphology, 
statistical tests, structural parameters, mass- function slope, 
and comparison with nearby OCs) that both objects are old 
star clusters. 

FSR 1716 is significantly absorbed (Ay ~ 6.3) and 
projected not far from the bulge. Its field-decontaminated 
CMD morphology is very similar to that of the ^ 7 Gyr 
well-known OC NGC 188. Indeed, its CMD can be well 
represented by isochrones with ages older than ^ 6 Gyr, 
both of Solar and sub-Solar {[Fe/H] ^ —0.5) metallic- 
ity. We adopted the former as the metallicity of FSR 1716, 
because of its relatively central location. Alternatively, we 
cannot rule out the possibility that FSR 1716 is a low-mass, 
loose (Palomar-like) globular cluster. FSR 1716 is located 
inside the Solar circle, « 0.6 kpc (in the case of an OC) or 
K. l.Skpc (GC). 



The CMD morphology of Cz 23 is consistent with that 
of the ~ 4 Gyr old OC M67. Similarly to FSR 1716, it 
can be well represented by Solar and sub-Solar metallic- 
ity isochrones, but with ages in the range 4 — 6 Gyr. With 
the w 5 Gyr and Solar metallicity solution, we find that 
Cz23 is projected nearly towards the anti-centre, located 
K, 2.5 kpc outside the Solar circle. 

The core and cluster radii of FSR 1716 and Cz 23 are 
small when compared to a set of open clusters in the Solar 
neighbourhood. However, such radii are comparable with 
those of other OCs of similar old age. Besides, the mass 
functions appear to be much flatter than Salpeter's IMF, 
especially FSR 1716, which seems to present an increasing 
depletion in the number of low-mass stars. As a conse- 
quence, the total masses presently stored in stars in both 
clusters are lower than ^ 200 Mq. Such low values prob- 
ably reflect the several Gyr-long period of mass loss due 
to stellar evolution, tidal interactions with the bulge (pos- 
sibly in the case of FSR 1716), disk and giant molecular 
clouds. Actually, because of its low mass content and flat 
mass function , Cz 23 may be evolv ing into an open cluster 
remnant (e.g. iPavani k Bicall200l . 

Comprehen sive catalogues of star cluster can didates, 
such as that of IFroebrich. Scholz k, Raftervl (|2007D . should 
be further explored with fleld-star decontamination algo- 
rithms and other tools, so that the nature of the candidates 
can be probed and the age derived. It is remarkable how 
the decontamination tool unveiled the intrinsic CMD se- 
quences of FSR 1716, separating it from the crowded field 
population. Consequently, the characterisation of FSR 1716 
and Cz 23 as OCs older than ~ 4 Gyr represents an impor- 
tant increase (« 10%) to the known population of such 
objects in the Galaxy. In particular, FSR 1716 is the most 
recent addition to the 8 open clusters older than 6 Gyr so 
far identified (WEBDA). In this sense, Cz23 (FSR 834), 
and especially FSR 1716, can be considered as rare exam- 
ples of extreme dynamical survivors in disk-regions where 
most open clusters are short-lived. 
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